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Abstract 

A generic class of string theories predicts the existence of hght moduh fields, 
and they are expected to have masses comparable to the gravitino mass 
777.3/2 which is in a range of 10~^keV-lGeV in gauge-mediated supersymmetry 
breaking theories. Such light fields with weak interactions suppressed by the 
Planck scale can not avoid some stringent cosmological constraints, that is, 
they suffer from 'cosmological moduli problems'. We show that all the grav- 
itino mass region 10~^keV ^ ~ IGeV is excluded by the constraints 
even if we incorporate a late-time mini-inflation (thermal inflation). How- 
ever, a modiflcation of the original thermal inflation model enables the region 
10~^keV ^ ~ 500keV to survive the constraints. It is also stressed that 
the moduli can be dark matter in our universe for the mass region 10~^keV 
<m^< lOOkeV. 
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I. INTRODUCTION 



Supergravity theories, which describe the low energy dynamics of string theories, gener- 
ically have a large number of fiat directions in their field spaces We call the scalar 
fields corresponding to these fiat directions moduli fields, or moduli, simply. Moduli are 
expected to take their values of the order of the Planck scale Mpi ~ 1.2 x lO^^GeV in 
the very early universe, because it is the only scale appearing in supergravity actions. At 
the later epoch in the universe's evolution supersymmetry (SUSY) is spontaneously broken. 
SUSY- breaking effects may lift the fiat potential for moduli through some non-perturbative 
dynamics and generate their masses comparable to the gravitino mass : ~ ""^3/2 0]- 

Moduli fields are generally accompanied by different kinds of cosmological problems 
depending on the values of their masses. These problems are divided into two classes dis- 
criminated by if or not the moduli decay takes place in the early universe. That is, if it does 
the radiation produced by the moduli decay may conflict with some cosmological observa- 
tions, and even if it does not a tremendous amount of moduli energy density itself causes 
disaster. Let us make a crude estimation of moduli lifetime and determine its dependence 
on m^. From dimensional analysis the decay width of the moduli is roughly 

Mil' ^' 
where N denotes the number of decay channels. The lifetime is given by 

r, = r-.lxlO-=eciV-'(^)-'. (2) 

In conventional hidden sector models, the gravitino masses m^i2 are about ITeV, so 
from eq.(|D we see that the moduli lifetime is much shorter than the age of the present 
universe ~ 3 x 10^'^sec for ~ Tn^i2- Therefore, in these models one must worry about 
whether or not the radiation produced by the moduli decay spoils the scenario of big bang 
nucleosynthesis. The reheating temperature derived from the width (|ip is 

3 

Tn ~ ^MpiT^ ^ O.SMeV iV^ (j^) ' • (3) 

In order for the big bang nucleosynthesis to work, high enough reheating temperature is 
needed (T^ ^ lOMeV), and this requirement constrains the moduli mass to be lower bounded 



m„ 



> 



0(10) TeV iQ. 



However, this relatively large moduli mass m^{r^ ^3/2) ~ lOTeV is realized only in a 
specific class of hidden sector models. Thus, one is faced with a difficulty that the process of 
nucleosynthesis does not proceed efficiently in a generic class of hidden sector models . 



^ It might appear that a large amount of entropy production by the moduh decay predicts an 
extremely smaher value of the baryon-to-entropy ratio than the one observed today even for ^ 
10 TeV, but the Affleck-Dine mechanism for baryogenesis M naturally explains the present value 
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Lyth and Stewart proposed a mechanism so-called thermal inflation |^ which solves this 
problem by diluting extensively the cosmic energy density of the moduli and consequently 
decreasing the number of photons produced by the moduli decay sufficiently so as not to 
upset the nucleosynthesis. 

Gauge-mediated SUSY-breaking models to which we pay our attention especially in 
this paper, have the gravitino mass in a range of 10~^ keV-lGeV. As can be seen from eq.(0), 
in the mass range 10~^keV^ < 100 MeV the moduli do not decay sufficiently fast and 
are left in the present universe. Generically, moduli fields are expected to take their vacuum- 
expectation values (vev's) of the order Mpi, which is so large that the moduli energy density 
easily exceeds the critical density of the present universe. Although the moduli lifetime (0) 
tells us that a class of these models with 100 MeV ^ ^ 1 GeV does not suffer from the 
moduli problem of their too much amount of energy density, a constraint from the cosmic 
7-ray backgrounds is crucial in this mass region [§,0 . 

Our main task in this paper is to apply the thermal inflation mechanism to gauge- 
mediated SUSY-breaking models 0, and argue the possibility that these models could pass 
the above cosmological constraints. We assume ~ rn3/2 throughout this paper. 

Our analysis consists of two parts. In the first part (sections 2 and 3) we adopt a thermal 
inflation model proposed by Lyth ane Stewart In this model there is a pseudo Nambu- 
Goldstone boson (called R-axion) arising from a spontaneous breakdown of an accidental R 
symmetry. We show that the mass of R-axion is always much smaller than the flaton mass 
in the gauge-mediated SUSY-breaking models. Thus, the flaton decays mainly into these 
R-axions and as a consequence the dilution of the moduli's energy density becomes much 
milder. We show that there is no parameter region surviving the cosmological constraints 
in this type of thermal inflation model. 

We modify the above thermal inflation model, in the second part (section 4), by adding a 
small explicit breaking term of the R symmetry to suppress the R-axion decay of the flaton. 
We find that a parameter region for the moduli mass, 10~^keV^ m<^ ^ SOOkeV |]TU[, survives 
the cosmological constraints. It should be noted here that the Affleck-Dine baryogenesis 



does work for these gravitino mass region [11|. We also stress that the moduli themselves 



could be the dark matter in the present universe for the region, 10~^keV^ ^ 100 keV. 
This may be a very crucial observation, since there is no dark matter candidate beside the 
moduli themselves because of the substantial dilution of any relic particle abundance by the 
late-time thermal inflation 0. The last section is devoted to discussion and conclusions. 



II. THERMAL INFLATION 

In this section we review a thermal inflation model which was proposed by Lyth and 
Stewart 0. Inevitable ingredients of the model are flaton fields, which are characterized by 



This attempt has been first proposed in Ref. [11|. 



Peccei-Quinn axion with high values of the decay constant Fa — 10 -10 GeV may be another 
candidate for the dark matter |12|. See also Ref. |jl 
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their 'almost flat' potentials and 'large' vacuum-expectation values (vev's). Here, 'almost 
flat' means that the roots of the second rank derivatives of potentials at their minima, 
namely the masses of corresponding flaton particles, are of the order of SUSY-breaking scale 
and the flaton flelds have 'large' vev's which are much larger than the scale. For simplicity, 
let us concentrate on the case in which only one flaton fleld exists, even though more than 
one flaton cases are also valuable to analyze. 

In order to guarantee the flatness of potential, we assume that the flaton's potential pos- 
sesses some exact (or at least approximate) global symmetries. If not, nonrenormalizable 
higher order interactions would induce a tadpole term through the SUSY-breaking effects. 
One familiar example of such global symmetries is U{l)ii symmetry, but it should be explic- 
itly broken to some discrete symmetry, say Z„, by a constant term C in the superpotential 
which is required to cancel the cosmological constant. Thus, in this paper we postulate that 
the superpotential for the flaton exhibits Z„+3 (n > 1) symmetry and takes a form 



W = C + E (^ + 3)^ Mi"+')''^'' 



where X is the flaton chiral superfleld and (Ai = 1) coupling constants. M* denotes the 
cut-off scale of the models we consideif]. As will be verifled later the vev of the superpotential 
(^) is dominated by the constant term C, and hence we take 

\C\ (5) 

to cancel the vacuum energy. Here, Mq is the reduced Planck scale, Mq = Mpi/\/Stx ~ 
2.4 X lO^^GeV. 

Then, the effective potential of the flaton X at low energy scale is represented as 

V.,AX) - V„ - ™S|Af + ^^(X"« + A—) + ^|A-|-«, (6) 

where we have used the same letter X for the flaton complex scalar fleld as for the corre- 
sponding superfleld, hoping that readers should not be confused. Note that the dynamics of 
the flaton fleld X is governed by the leading term in the superpotential (1) for X -C and 
we have neglected higher order terms in eq.(^ since the vev of the flaton is much smaller 
than the cut-off scale M*. Vq is determined by the requirement that the cosmological con- 
stant vanishes Pi The quadratic term in eq.(p|) is induced by SUSY-breaking effects WA 



^ Even if one assumes an approximate U{1)r symmetry instead of Z„4_3, one reaches the same 
conclusion as in this paper, since the higher power terms (/c = 2, • • • , oo) in the superpotential 
are practically negligible as seen below. 

^ In the original paper Q the cut-off scale is taken to be at the Planck scale. Here, we regard 
it as a free parameter to make a general analysis. 

^Vb is at the tree- level given by 
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and we assume the mass squared at the origin X ~ to be negative and of the order of 
SUSY-breaking scale. The vev of the flaton is given by 

(X) ^ M ^ (-l^)'™(moM:)^, (8) 
where we have neglected 0(^3/2/^^0) terms. (Remember that the formulae given below are 



also vahd only up to 0(7/13/2/^0). See the paragraph following eq. ([T8|) .) Then Vq is given 
by 

V, ~ HlmlM^ (9) 

" n+2 ° ^ ' 

Here, reflecting the Z„+3 symmetry the degenerate minima cause a potential domain wall 
problem. We will come back this point in section 4. 

The superpotential (^ has an approximate U{1)r symmetry if we neglect the higher 
order terms and then the imaginary part of X becomes a pseudo-Goldstone boson called 
R-axion. The R-axion receives a mass from the constant term C which breaks the U{1)r 



symmetry explicitly [15]. If we parameterize X as 



X^(-L,^M)exp(-|^). (10) 

with X representing the real flaton field and a the R-axion, then the flaton and the R-axion 
mass squared are estimated a^ 

mj~2(ra + l)m^, (11) 
2 n{n + 2>) n(n + 3) 

+ 2 ^2{n + l){n + 2) 

The thermal inflation occurs if the flaton field does not sit at the true minimum of the 
potential but at the origin in the early universe. To realize this initial condition for the 
thermal inflation the flaton must have interactions with other fields in the thermal bath 
of the universe^. From the finite temperature effects, the effective potential in the early 
universe takes a form as 



Vo-|F|^--^|C|^ (7) 
where F is the SUSY-breaking F term. 

^ Here we choose the phase of the constant term in eq.(^ as C ~ so that the R-axion 

mass squared is positive. 

^ For example, a Yukawa interaction as VF = gX^^ is sufficient. The fields ^ and ^ receive a mass 
~ gM when the flaton sits at the true minimum, but for the flaton field values near the origin 
the fields ^ and ^ become light and could be in the thermal bath if they couple to the particles in 
the standard model and the temperature T is larger than the mass m^. 
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where c is a constant of 0{1) and T is the temperature of the universe. Then, at high 
temperature T ^ Tc ^ mo the flaton field sits near the origin and produces the vacuum 
energy Vq. 

If the energy of the universe is dominated by the radiation just before the thermal 
inflation, the vacuum energy Vq becomes comparable to the radiation energy at the cosmic 
temperature ~ K)^''^- Thus, for Tc ^ T ^ T^,, the vacuum energy of the flaton field 
dominates the energy of the universe and a mini-inflation, i.e. thermal inflation takes place 
0. On the other hand, if the moduli oscillations dominate the energy of the universe before 
the thermal inflation, the temperature at the beginning of thermal inflation is estimated as 

~ {Vq / {m^Mc)Y^^ . Here are the moduli masses. 

In considering the history of the universe after the thermal inflation, the flaton decay is 
crucial since it is most responsible for the entropy production. Here we shall list the relevant 
decay channels and compute quantitatively the decay rate for each channel which is needed 
to trace the physics following the thermal inflation epoch. 

The only possible renormalizable interaction of the flaton with SUSY standard model 
particles is 

W = XXHH, (14) 

where H and H are Higgs chiral supermultiplets. When the flaton field develops the vev, 
{X) = M, the Higgs multiplets acquire a mass AM, and for this to be at most the electroweak 
scale the coupling constant A has to be set so small, A ^ fin/M. Here fin is the SUSY- 
invariant mass for the Higgs multiplets. If the flaton mass is large enough, the flaton decays 
into a pair of Higgs fields^. The decay width is represented as 

1 TTl'^ 

- '^'IfaMt. (15) 

where Ch is a constant parameter satisfying Ch ^ and we have ignored the masses 

of the Higgs fields. 

Even if the above decay process would be kinematically forbidden, the flaton could decay 
into two photons through one-loop diagrams |T^Q. The width of this decay channel is given 

by 



^ Here, Higgs field denotes a Higgs boson or a Higgsino. 

For example, the coupling Ch is Ch = (XM/m^)^ in the case the flaton decays into two Higgs 
bosons. And Ch = {XM/m^)'^ for the case the flaton decays into two Higgsinos. Thus in this 
analysis we shall assume that the coupling Ch is a free parameter of C/i < 0(1), since A ^ fj-n/M. 

This decay process is induced by the Higgs fields loop diagram through the interaction ( |l^ or 
by charged ^ loop through the Yukawa interaction which is required in order that the flaton sits 
around the origin during the thermal inflation. If ^ has a color charge, the decay into two gluons 
occurs and becomes a dominant decay for the flaton mass > a few GeV. 
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r , (16) 

Through the above decay processes, the flaton energy is transferred to the radiation and 
reheats the universe 0- 

One should note that the flaton can decays into two R-axions, if kinematically allowed. 
The decay rate is calculated as 

1 m?, , ^ 

Here, we have neglected the R-axion mass. The R-axions produced by this process succes- 
sively decay into two photons similarly to the flaton decay and its decay rate is estimated 

as 

In gauge- mediated SUSY breaking theories the predicted gravitino mass range (^3/2 — 
10~^keV-lGeV) indicates that the flaton decay into two R-axions is most likely allowed, 
since the flaton may obtain a mass of the order of the SUSY-breaking scale {m^ ~ lOGeV- 
ITeV) and the R-axion has a mass nia ~ ^/^n^jWh^ (^^^ 6qs.(|ll|) and ([T^) ) Thus, in the 
following, we consider the case that the flaton can decay into standard model particles (Higgs 
flelds or photons) and two R-axions. (For the case that the R-axion decay is forbidden, see 
section 4.) 

Now we are ready to describe the history after the thermal inflation. Our aim here 
is to estimate the entropy produced by the flaton decay processes. At the temperature 
T ~ Tc the flaton fleld starts to roll down to the true minimum of the potential (^) and then 
oscillates around it. When the Hubble parameter becomes of the order of total width of the 
flaton, the flaton x decays into both standard model particles and R-axions. The energy 
of R-axions can not be transferred to the radiation at this time because it has only weak 
interaction suppressed by 1/M with particles in the thermal bath. Then only the energy of 
the standard model particles is transfered to the radiation to reheat the temperature of the 
universe Tji^sm at the flaton decay epoch. The ratio of the entropy densities just before to 
after the flaton decay is estimated as 

Asm ^ l + {l-ea)^ ^ ^° , (19) 

where g^^ is the effective number of degrees of freedom and denotes the branching ratio of 
the flaton decay into two R-axions. 



A possible decay into two pairs of bottom and antibottom quarks is strongly suppressed by 
phase volume effects. 

In Ref. hidden sector models of the SUSY breaking are considered where the gravitino has 
a mass 7713^2 — w-^ — 100 GeV-1 TeV. In this case the R-axion has a mass of the order (see 
eq.(|l2|)) and the flaton decay into R-axions may not be allowed. 
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Just after the flaton decay, the energy density of the R-axion per the entropy density is 
(27rV45)g*r^ ) a1^' '^^^ energy densities of the R-axion and radiation are both diluted by 
the expansion at the same rate where R is the scale factor of the universe. But after 
the R-axions become non-relativistic particles, their energy density pa is diluted at and 
dominates over the energy of the radiation unless ~ 0. 

The R-axion decay into two photons occurs at the Hubble parameter H ~ ^a^2-y and 
the universe is reheated again with the temperature Tr. At this time, the R-axion decay 
increases the entropy by a factor 

4: Vo 1 (2ma\ 

^'^ ^ ^ + '"3(27rV45)^.T3rKA5M J ' ^ ^ 

Then, all the vacuum energy of the thermal infiaton (i.e. fiaton) is eventually transferred 
to the radiation. The reheating temperature Tr is estimated from eq. (|l^) as 

Tn ^ 1.1 X 10-^ ^'J'M'f 

And the thermal inflation increases the entropy of the universe by a factor 

A - . A. ^ 1 + (1 - e.)3(2,./45)^^^3y^^^^^ + ^'^g (2^^Vi%jf^ l^j ' ^ ^ 

As we have described the thermal inflation reproduces the entropy at the late-time of 
the universe's evolution which dilutes substantially the cosmic energy densities of any relic 
particles such as moduli fields. The dilution factor is given by eq. (^2]) when the two R-axion 
decay of the flaton is allowed. The dilution factor in the other case will be given in section 
4. We should note that the requirement Tr >10 MeV in eq.(pT]) leads to the upper bound 
for M which justifies our ansatz in eq.(||). 



III. MODULI PROBLEM IN GAUGE-MEDIATED SUSY BREAKING THEORIES 

In this section, we consider the cosmological moduli problem in gauge-mediated SUSY 
breaking theories. The predicted gravitino mass m3/2 — 10~^ keV-1 GeV indicates that the 
string moduli fields may have lifetimes longer than the age of the universe. Then one clear 
problem arises that the energy densities of the moduli overdose the universe. Moreover, it 
has been pointed out recently that the contributions to the cosmic X(7)-ray background 
from the moduli decays are very dangerous P,p!0|. The only known possibility to solve these 



problems is the thermal inflation. Therefore, we examine whether the thermal inflation 
discussed in the previous section can solve these moduli problems. 



A. The energy density of string moduli 

First, we discuss the problem of the cosmic energy density of the moduli fields. To derive 
conservative cosmological constraints, we assume that at least one modulus field exists with 
a mass m<^ ~ m3/2 and it has an initial value of the order of the gravity scale Mq- (Here, we 
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have chosen vev of the modulus field at the origin). Generalization to the case that many 
moduli fields exist is straightforward. 

Let us show by exphcit calculation that the modulus energy density exceeds substantially 
the critical density of the present universe if the modulus is stable. When the expansion 
rate of the universe becomes of the order of the the modulus mass, H ~ m^, the modulus 
field (f) starts to oscillate around the minimum of the potential with the initial amplitude (^o 
of the order Mq- At that time the energy density of the modulus coherent oscillation is 

P<t> = \f^l4>l. (23) 

and the energy density of radiation is of the same order of p^. Thus, the temperature of the 
universe when the modulus starts to oscillate is estimated as 

/ 90 



Mcm^, (24) 

.T.2xlO»Gev(^)'^ (25) 

Here, we have assumed that the modulus coherent oscillation begins after reheating of the 
ordinary infiation completes, and that its reheating temperature is higher than T^]^ Then 
the energy density per the entropy density is given by 



s (27rV45)^7*T|' 



.O.x.O»GeV(^)-(A^. (27) 

This ratio takes a constant value until the present if no entropy is reproduced, because the 
densities of the energy and of the entropy are diluted at the same rate as as the scale 
factor R increases. 

On the other hand, the critical density of the present universe is given by 

^ = 3.6 X 10-9 h'^ GeV, (28) 
s 

where h is the present Hubble parameter in units of 100 km/sec/Mpc. When we consider 
gauge mediated SUSY-breaking theories where the predicted modulus mass is ~ "^3/2 — 
10~^ keV-1 GeV, one can see from eq. (^7[) that the energy density of the modulus coherent 
oscillation overdoses the universe if the modulus is stable until now. The thermal infiation 
increases the entropy of the universe by the factor A as shown in eq.(p2D and dilutes the 
energy density of the modulus given by eg. ([2B|) . 

Now we adopt the thermal infiation model in the previous section and estimate the 
minimum value of the present energy density of the modulus field 0. The relevant dynamics 



2 



14 



In some models of inflation, such as a chaotic or hybrid inflation, one may easily have the 



reheating temperature higher than in eq.(25) |16| 
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is determined by two mass scales mo, and the branching ratio ea for a given gravitino 
mass 1713/2 (— rriif,). Since A in eq. (|2^) takes its maximum at = 1 as far as uia -C m^, 
we put = 1 to obtain the most efficient dilution factor In the following analysis we 
take two free parameters and Tr instead of mo and M^, and search the minimum energy 
density of the modulus in the present universe. 

The amount of the present energy density of modulus takes different forms depending 
on whether the modulus field begins to oscillate before the thermal infiation or after the 
end of that. For the moment we assume that the modulus coherent oscillation starts before 
the thermal infiation, i.e. m^ > Hti- Here, Hti — ^/Vo/ {\/3Mg) is the Hubble parameter 
during the thermal infiation. Then, from eq.(p6|) the present energy density of such a 
modulus ("big-bang" modulus) is 

PA ^ ^^y^ 1 (29) 



s Jbb (27rV45)(?,r| A' 

Furthermore, the modulus energy is reproduced after the thermal infiation. Because, during 
the thermal infiation, the modulus sits at the minimum of the potential which is shifted from 
its true vacuum by an amount of 50 ~ (Vo/m^M||)0o and restarts to oscillate around 
the true minimum with an amplitude 6(f) after the end of the thermal infiation. Then the 
present energy density of this "thermal-infiation" modulus is estimated as 

•nim-/2 1 ,3„, 



s Jti (27rV45)(7,T3 A' 



Therefore, the present total energy density of the modulus with mass m^ > Ht.i. is iToH 



max 



S / BB \ S / TI 




(n + !)(. + 2)^ (nV (<k_y_rl_ 



^(n^3V I m_ / \ I /i^l/2 3/2 

G '"3/2 



Here, notice that it depends only on Tr and not on m^ since Tc ~ m^. The reheating 
temperature should satisfy Tr ^ 10 MeV to maintain the success of the big bang neucle- 
osynthesis. Then if we choose the minimum value of Tr ~ 10 MeV in eg. (PTf) , we obtain the 
lowest value of the present modulus energy density for ~ m^ and 0o — Mq, 

On the other hand, if the modulus mass is smaller than Hti, then the modulus oscillation 
begins after the end of the thermal infiation with the amplitude 0o ~ Mg- The entropy of 
the universe when the modulus starts to oscillate is estimated as 



° When the flaton decay into Higgs fields is forbidden, we have 6^ — 1. When it is allowd, 6^ 
depends in eq.([l5|). However, we find that the dilution factor A takes the maximum value at 
Ch — which means — 1. 



10 



where we have used the fact that the entropy is diluted by the expansion of the universe at 
the rate . Then the present abundance of the modulus is[3 



{- 



15 \ 02^0 1 



_ 1 [2(n + l)(n + 2)]^/" m^^TR 
^16 Inin + 3) U'^cj mj/' 



(34) 
(35) 



''3/2 

In this case, the mass of flaton x, (since m,^ < i^T/) should satisfy 



(n + 2)V^(n + l)3A4 
> 3.5 X 10 ^^^^ ^ 3^j3/7 7713^3 Tj{ , (36) 

and this gives the lower bound as 

^-j^^l.6GeV [^(^^3)p/, [-^j , (37) 

for the minimum value of the reheating temperature Tr ^ 10 MeV and 00 — Mfj. Here we 
have neglected the region where the gravitino mass is less than 

2 ji8/5 

m3/2 ~ 2.4 X 10^ 7 TTTTTTr^ TTTTT^— Itf, (38) 

^' (n + 1)1/10 [ri(n + 3)]6/5]VfJ/5' ^ ^ 

since the vev of the flaton M exceeds the cutoff scale M^, and our effective treatment of the 
flaton potential breaks down. 

To compare our result with the critical density of the present universe, we show the 
obtained lower limit for 

Pc 

in Fig.l. In this analysis we have assumed that the modulus is stable and hence this figure 
shows the cosmic energy density of the stable modulus' coherent oscillation. We find that 
it exceeds largely the critical density of the present universe. Notice that in this figure we 
take the case of n = 1. 

So far, we have assumed that the modulus field is stable. However, it is not valid, since 
the modulus field may couple to the ordinary particles through some nonrenormalizable 
interactions. The most plausible candidate for the modulus is the dilaton in string theories 



1^ Here the modulus oscillation is assumed to start at least before the R-axion decays, i.e., 
Va-.2-y ^ 1.1 T|/Mg ~ 4.5 X 10-23 GcV {Tr/W MeV)2. 
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and it decays most likely into two photons: — > 27 [|] through a coupling to two photons 
written as 

= («) 

Here, we have introduced a dimension-less parameter h which depends on the type of super- 
string theories and compactificationsQ. In the present analysis, we take 6 as a free parameter 
of the order one representing the various compactifications in string theories. Then, the life- 
time of the modulus is given by 

~ — f ~ 7.6 X 10^3 sec - . (41 

Thus if ^100 MeV the modulus would decay within the age of the present universe 
and its energy would be diluted below the critical density of the universe. But for such a 
case another stringent constraint should be considered. 



B. Constraint from the cosmic X(7)-ray backgrounds 



Even if the lifetime of the modulus is longer than the age of the universe, the modulus 
particle decays into photons in the past universe. Thus, the produced radiation will con- 
tribute to the cosmic X(7)-ray backgrounds and the observed backgrounds give a constraint 
on the mass and the lifetime of modulus [pHlO|. 

The photon number flux induced by the modulus decay is given by M 



dt -^{l + z)5 {E^il + z)- 




exp 



3/2 



(42) 
(43) 



where to is the age of the universe, z the red-shift. Ho the present Hubble parameter and 
the energy of X(7)-ray. And n^pfl denotes the present number density of the modulus if it 
would be stable. Here, we have assumed that the present total density parameter is f2o — 1. 
The detailed derivation of the photon flux in eq. (|42|) will be given in Appendix. 

We may obtain a constraint on Q^h"^ [Q^ = m^n^^o/ Pc) by requiring that the maximum 
value of the flux in eq.(^) should not exceed the observed X(7)-ray backgrounds ||T9H2T 
The observational data are fitted by the following power-low spectra [H 



The modulus decay into two neutrinos is suppressed since it has a chirality flip and vanishes for 
massless neutrinos Similarly to the two photon decay the modulus may decay into two gluons. 
For such a case the lifetime of the modulus may be shorter than that in eq. (41) by a factor about 
9. 

For example, the dilaton has a coupling b = \/2 for a compactification of the M-theory 
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7,of)s 



{E. 



7/ 



cm -sr-sec 



8 (E^/keV)-o-4 
380 (E^/keV)-^-^ 
2 (E^/keV)-o-^ 
1.6 X 10-2 {E^/MeVy 
1.5 X 10-3 {E^/MeVy 



-1.8 



O.lkeV < < 25keV 
25keV <E^< 350keV 

350keV <E^< IMeV 
IMeV < < 20MeV 

20MeV <E^< lOGeV 



(44) 



< 

r\-/ 



10 GeV 



The result is also shown in Fig.l. We find that all mass region 10"^ keV ^ 
is excluded completely by the observed X(7)-ray backgrounds. 

In summary, we have shown that the thermal infiation could not dilute the energy density 
of the modulus sufficiently so as to lower it below the critical density of the present universe, 
if the modulus is stable. Moreover, the observed cosmic X(7)-ray backgrounds put a more 
stringent bound on the modulus with mass 200 keV ^ < 10 GeV even if modulus is 
unstable. Therefore, the cosmological string moduli problem is not solved by the thermal 
inffation for the moduli masses less than 0(10) GeV. This excludes the mass range predicted 
in gauge-mediated SUSY-breaking theories as long as ~ 



IV. THERMAL INFLATION WITHOUT THE R-AXION DECAY AND THE 
COSMOLOGICAL MODULI PROBLEM 

A. Modified thermal infiation model 

If one sees carefully our scenario of the thermal history described in the previous section, 
one might recognize that the reason why the thermal inffation mechanism has failed to 
solve the moduli problem could be attributed to the R-axions produced by the ffaton decay 
1^. Indeed, the energy density of the relativistic R-axions decreases faster than that of non- 
relativistic particles, and as a consequence the R-axion decay into photons releases much less 
entropy than in the case that the decay mode of the ffaton into R-axions is not open. Thus 
we might expect that some mass regions in gauge-mediated SUSY breaking theories could 
survive the cosmological constraints if we could improve the model to forbid energetically 
the ffaton decay into R-axions. 

Besides the stringent constraints discussed in the previous section, there is another dif- 
ficulty, i.e. the domain wall problem, which we have ignored so far. The origin of this 
problem is the degenerate minima of the potential which possesses the discrete symmetry 
Zn+3- Thus, in order to eliminate the domain walls we have to add a small term to the 
potential which breaks the discrete symmetry explicitly. 

One of economical modifications of the model which satisfies both of the above two 
requirements is, for instance, to add a linear term in the superpotential 

6W = aX, (45) 
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In hidden sector models for the SUSY breaking the flaton decay into R-axions may not be 
allowed energetically (see footnote 13). Even if it is allowed the R-axions decay just after the 
flaton decay, since the R-axion mass is of the order m^, which reheats the universe immediately. 
In any case sufficiently large entropy is produced and the moduli problem may be solved 



13 



which breaks the Z^+s symmetry completely together with the constant term C. To collapse 
the domain walls before its energy dominates the universe, the dimensionful parameter a is 
required to be ET 



\a\ 



(46) 



The above explicit breaking term (^5|) in the superpotential modifies the low energy 
potential of the fiaton in eq.(BD as 



aC 

K// (X) = Vo-2j^{X + X*) - ml\X\ 



G 



a 

Mr 



*n+2\ 



n 



C 



n + 3 M^Ml 



*n+3N 



+ 



IX 



2n+4 



(47) 



In the following, instead of a we use a dimensionless parameter x defined by a = 
— xiVf^+^/M^, for simplifying the expressions below Then the vev of the fiaton is given 
by 



(X) = M 



1 



(n + 2)(l -x) 



2{n + l) 



(moM: 



n + l 



and Vn is 



nil — X) + 1 n , ^9 



(49) 



(50) 



It should be noted that the explicit breaking term (^Sf) does not affect much the dynamics 
of thermal inflation. In the early universe the potential of the flaton near the origin takes 
the form as (see eq.(|T3|).) 



V;//(X)^(cT2-m, 



X 



2am3/2 



cT^ — ml 



cl — mi 



(51) 



Although at the high temperature T ^ Tc ^ mo the flaton does not sit at the origin, the 
position of the minimum is 



(X) 



2am3/2 
cT^ — m, 



2x 



■0 ^(n + 2){l-x) 



cT'^ — mn 



M, 



(52) 



We can write down the relation of x and a in terms of the parameters appearing in the 
superpotential as 



a 



, n + l Jl/fn + l 



[(n + 2)(l 

up to the terms of O (771,3/2 /?7io)- 



^mS^'M,-+\ (48) 
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and the deviation from the origin is suppressed by a factor of C(m3/2/'mo) compared to the 
vev M for the true minimum 0, as far as cT^ — ^ m^. Furthermore, the deviation of 
the vacuum energy of the flaton from Vq is estimated as 

Ax^ mo/9 

^ - 2 ^0, (53) 

which has a suppression factor of 0(171^2/ ''t^o) and is consequently neghgible. Therefore, we 
find that if we take x to be not so close to unity the dynamics of the thermal inflation is not 
modified much by adding the explicit breaking term p5|). 

Since the linear term (^) in the superpotential breaks also an approximate U{1)r sym- 
metry, the term gives additional contributions not only to the fiaton mass but also to the 
i?-axion mass. Then their masses are given by 



2in + l)-nx 2 /..N 
-mo, (54) 



X 



m„ ~ -YZT^^o- (55) 



Therefore, for the region^^ 



Xr 



< X « 1), (56) 
5n + 8 ^ ^' ^ ^ 



the fiaton decay into two i?-axions is kinematically forbidden. In the following we consider 
the thermal inflation model with the explicit breaking term (|45| ) satisfying eg . (^6|) , and re- 
examine whether the model can solve the cosmological moduli problem in gauge-mediated 
SUSY breaking theories. 



B. Cosmological moduli problem with the modified thermal inflation 

First we argue the thermal history after the thermal inflation. At the end of the thermal 
inflation (T ~ T^) the fiaton field begins to oscillate around the true minimum of the 
potential (|47|), and when the Hubble parameter becomes comparable to the total width of 



If we take the value of x to be very close to unity (x < 1 is always satisfied by the definition of 
x), the position of the minimum at the high temperature is far from the origin. Thus we consider 
the model with x not so close to unity. 

Although this seems to be a highly restricted region, one can see by taking the original 
parametrization by a that it actually corresponds to a broad one, namely 

, , 2 n + l / 5n + 8 \^(^ ..7^ ^ 
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the flaton, the flaton decays only into SUSY standard model particles since the decay into 
i?-axions is not allowed. As discussed in section 2, the flaton decays dominantly into Higgs 
particles if kinematically allowed, or when the flaton mass is smaller than the threshold of the 
decay into two Higgs particles the flaton decays into two photons. Each width is represented 
in eqs.(p^) and (|1^) , respectively. In both cases, the flaton energy is transferred to radiation 



and reheats the universe immediately. Then the modifled thermal inflation model increases 
the entropy by a factor given by putting ea = and denoting Tr^sm by Tr in eq.(|T^, 

^"^ + ^(27r745^lT3T«' ^^^^ 

where the reheating temperature Tr is obtained by the decay width of the flaton. 

For the case that the flaton can decay into Higgs particles (m^ > 130 GeV) 0, Tr is 
represented as 

T^^0.14F(C,)^^^, (58) 



where F{Ch) is deflned by 



"em 



477 



1/2 

(60) 



Thus in this pointed out in Ref. [0, the reheating temperature Tr can be taken to 

be high enough (Tr > 10 MeV) to maintain the success of big bang neucleosynthesis. On 
the other hand for the case < 130 GeV, the reheating temperature is estimated from 
eq.(|TB|) as 

Tr ~ 1.1 X 10-^ ^^"^ , (61) 

" (loo^v) [-^)- ^''^ 

Now let us estimate the cosmic energy density of coherent modulus oscillation in the 
present universe. The present energy density of the "big-bang" modulus is given by replacing 
A in eq.(P^) by the expression of the entropy production (0), hence 

, , , , — -. (63) 




23 Prom the experimental lower bound on the mass for Higgs bosons |23|, we take the Higgs mass 
to be 65 GeV in order to obtain as conservative constraints as possible. 
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On the other hand, the present energy density of the "thermal-inflation" modulus is esti- 
mated from eq. (|30|) as 

^) -0-38 f#V4^. (64) 

sJti \MgJ mlMl ^ ' 

Then we turn to estimating the lower bound of the total energy density of the modulus. 
The lower bound is given by 



~ max 
s y 



S / BB \ s / TI 



> 



—) 

S J BB \ S ) Tl' 



1-2^ ^ (65) 



Here, the equality is satisfled when {p^/s)bb = {p<i>/s)ti, i-e. when 

\ 3/''' 

' ^ M'J'^Tf, (66) 



2.5 X 10^ / Tc ^ 



~ :r-pi — I I m. 



where Cy^ = (^n+2)^{n+i)-nx] have used eq. (|6l|) for an expression of Tr, assuming 

< 130 GeV. Then, the lower bound of the moduli energy is estimated from eq. (|65D as 

PA ^ 4.8xio3 /0oy/T.y^/^ r^^/^ 

s)o^ \Ma) \mj mf'uT' ^ ^ 

Thus the lowest reheating temperature Tr ~ 10 MeV and x = Xmin give the minimum 
abundance as 



flrhh'^ > 3.9 X 10^^ 



3n + 8 



(68) 



for 00 — Mg and Tc ~ m^. However, if we put Tr ~ 10 MeV in eq.(pB|) we observe that the 
assumption < 130 GeV holds only when 

m^<m,^ 5.2 x 10"^ GeV (^)'^' . (69) 

Thus, for ^ rric the theoretical lower bound (|68D is realized in the case < 130 GeV 
and the main decay mode is x 27. 

On the other hand, for the modulus mass > rric we can obtain the minimum abun- 
dance lower than the r.h.s. of eq.(|H5|) by making use of eq.(p5|) which is applicable for > 
130GeV instead of eq.(^) 0. Therefore, for > rric the condition that the lower bound 
is saturated in eq.(B3) is 



In this case 777-3/2 > rric. Then, the flaton decay into Higssinos should be forbidden, otherwise 
the gravitinos produced in succesive decay of Higgsinos overdose the present universe. 
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4.3 I Tc \ ^ 5/14,^1/14^^4/7 

The lowest possible values for the flaton mass and the reheating temperature satisfying 
eq.(|70D are practically given by ~ 130 GeV and Tr = 10 MeV, yielding the minimum 
abundance in eg .(1651) as 



s/o \Mg) \m^J \ ) 

Comparing with the present critical density in eg. (]28|) , 0o — and ~ leads to 



3/4 

17^/1^ > 8.1x10-5(^^1 . (72) 

Tflfj) 



I GeVy 



In Fig. 2 we show the lower bound for the enegy density of the modulus predicted from 
egs.(^) and (0). A remarkable conseguence which distinguishes itself from the result in 
section 3 is that for all the moduli mass region 10~^ keV < < 10 GeV, the predicted 
lower bound can be taken to be below the critical density of the present universe. 

So far, we have been considering the case > Hti- Then let us discuss briefly what 



happens in the case < Hti- From egs.(0) and (|57|), the present abundance of such a 
modulus is given by 



1 



8 \Mg 



Tr, (73) 



for any fiaton mass. Here we have assumed that the modulus mass is larger than the 
decay width of the fiaton (See the footnote |T6[). Comparing with the case > Hti, the 



abundance in eg.(^) is always greater than the minimum abundances in egs. (|68[) and (|72D 
for ~ 10-2 keV-1 GeV. 

In Fig. 2 we also show the constraint from the observed X(7)-ray backgrounds which is 
derived in section 3. Then we can see that it excludes the modulus mass region 500 keV 
^ 1TL^ ^ 10 GeV. Therefore, we conclude that the cosmological problems are resolved in a 
class of gauge- mediated SUSY-breaking models with the gravitino mass 10-^ keV ^ m3/2 ^ 
500 keV, provided that the modified thermal inflation takes place and the fiaton decay into 
R-axions is forbidden. We should comment here that the ansatz in eg. (|^) is always satisfied 
for the parameter region we have analyzed. 



V. DISCUSSION AND CONCLUSIONS 

We have found in the first part of our analysis that due to the presence of too light 
R-axions the original thermal inflation model |^ could not resolve infamous cosmological 
moduli problems in gauge-mediated SUSY-breaking theories, and all the gravitino mass 
region, 10"^ keV ^ WI3/2 ^ 10 GeV, is excluded as long as ~ ^3/2 is fulfilled. 

In the next step, we have modified the original thermal inflation model so that the R- 
axion decay of fiaton is kinematically suppressed, and as a result we have succeeded to make 
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a region 10~^ keV ^ "^3/2 ~ 500 keV to survive the cosmological constraints. Moreover, it 
can be shown that a small window 1 GeV ^ "^3/2 ~ 10 GeV appears, if we take (po/Mo — 
0.01-0.1 that is nevertheless a slightly improbable condition]^. One should note, however, 
that this window is located in the region 7713/2 ^ MeV, in which it is far from easy 

to construct mechanisms to produce a sufficient number of baryons consistent with today's 



observational data |]TT],3 if one adopts the original Affleck-Dine mechanism [^]. Here, it 
may be interesting that an enough number of baryons will be created even for 7713/2 ^ 0{1) 
MeV if we adopt a variant type of Affleck-Dine baryogenesis discussed in Ref. p4[, since 



gauge-mediated models with 7^3/2 ~ 1 GeV are relatively easily constructed P5| , PB | 



In the present analysis we have adopted the effective potential (|g) for the flaton. We 
consider that the soft mass ml in eq.(|) is induced by the SUSY-breaking effects and hence 
the potential is applicable only for the scale below masses of messenger multiplets. 
Therefore, (X) = M should be less than 7r7^(~ 10^''~^^7r73/2 We have checked that 

this constraint is indeed satisfied for the mass region 10~^ keV ^ 7/13/2 ^ 10 GeV we have 
studied. (Note that (X) ~ 10^ GeV - lO^^ GeV for 7713/2 ~ 10"^ keV - 10 GeV.) 

We would like to emphasize that the moduli could be the dark matter in the present 
universe with their masses being possibly in the region 10~^ keV ^ ^ 100 keV. Fig. 3 
exhibits the region of the parameters and in which the conditions < pc and Tr ^ 
10 MeV are satisfied for = 100 keV. By noting that this region contains plausible values 
777^ ~ 10^-10^ GeV and ~ 10^^ GeV fortunately, one would be encouraged to expect 
that the above emphasis is by no means a nonsensical suggestion. The detailed study on 
the spectrum of the X(7)-ray emitted from the cosmic moduli will be given in Ref. p7[ . 



As concerns the moduli masses, it is another intriguing possibility that so-called "MeV- 



bump", an excess in the 7-ray background spectrum around ~ 1 MeV |20], might be a signal 



for the moduli decay M. Anyway, future analyses of X(7)-ray backgrounds in these mass 



regions are expected to provide us with more detailed information about the moduli . 

If the existence of light moduli fields is a generic prediction of string theories, it will 
probably present us with a plenty of subjects in low energy physics, typical examples of 
which are nothing but the cosmological moduli problems we have attempted to solve in this 
paper. 
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APPENDIX A: PHOTON FLUX 

Since the velocity dispersion of the moduli is negligible, two monochromatic photons with 
energy 7r70/2 are produced in the moduli decay. Thus, the spectrum S{E') of the photon 
per decay is written as 

S{E') = 2S{E' -m^/2). (Al) 
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This window will also appear if one adopts very small values for b ~ 0.01-0.1 in eq.(|4 
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The number density n,j){z) of the moduh at red-shift z is given by 

n^{z) = n^,o(l + zf exp{-t/r^). 
Then the present flux from the moduh is estimated as 

= ^ r dt'-n^,oexp{-t'/T^){l + z)26iE^{l + z) 
47r JO 



(A2) 



(A3) 



where to is the present time, the present energy of the photon and the factor {1 + z)~^ 
in the first equation represents the dilution due to the cosmic expansion. Here we have 
taken into account that the photon with energy had energy E' = {1 + z)E^ when it was 
produced at the decay time. The red-shift z is related to the cosmic time t by 



dt/dz = -H^\l + z)-'/^[no + (1 - fio - ^a)/(1 + + Qa/{1 + z] 



31-1/2 



(A4) 



where Hq is the present Hubble parameter, VLq the present (total) density parameter and 
the density parameter of the cosmo logical constant. Then photon flux is given by 



F.y(^E^ 



3/2 



rin^/2E^ 



X exp 



f{m^/2E^ 
1 



d{l + z)--—{l + z)-'^'f{l + z) 



where 



fix) = [fio + (1 - 1^0 - ^a)/x + nA/x^]-'/^. 

For ^0 = 1 and = 0, eq. (|A5|) is simplified as 



F^{E^) 



2ttt^Ho J 



3/2 



exp 



(2E^\ 



3/2 



) 



(A5) 



(A6) 



(A7) 



The flux F^ takes the maximum value Fr^^rnax at 



Er, 



2 

^ 3r^go ^2/3 



for > ^Hq ^ 
for < ^Hq'^ 



(A8) 
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FIG. 1. The lower bounds on the modulus densities O^/i^ in the original thermal inflation 
model for n = 1. The thick (thin) solid line represents the lower bound for the modulus density 
with > Hti {m^ < Hti)- The upper bound from the present critical density is represented by 
the dashed lines with a kink at ~ 100 McV. The kink appears because the constraint becomes 
weaker for ^100 McV due to the modulus decay. The experimental upper bound from the 
cosmic X(7)-ray backgrounds {b = 1) is shown by dot-dashed line. 
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n= 1 




m, [GeV] 



FIG. 2. The lower bounds on the modulus densities f2^/i^ in the modified thermal inflation 
model for n = 1 where the flaton to two R-axions decay is forbidden. The thick solid line with 
a kink at ~ 10 MeV represents the lower bound for the modulus abundance. The thin solid 
line for ^ 10 MeV represents the lower bound for the case the flaton to Higgs fields decay 
is forbidden. The upper bound from the present critical density is represented by the dashed 
line. The experimental upper bound from the cosmic X(7)-ray backgrounds (6 = 1) is shown by 
dot-dashed line. 
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FIG. 3. The contours of n<^/i^ in the modified thermal inflation model for n = 1 where the 
flaton to two R-axions decay is forbidden. We take the modulus mass = 10"^ GeV. The 
dashed line denotes Q^ph'^ = 0.25, and the dot-dashed line denotes = 1. In the shadow region 

Tr < 10 MeV in which the big bang nucleosynthesis is not operating well. 
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